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It is shown that optical synthesis of terahertz and millimeter-wave frequencies can be achieved us- 
ing two-mode and mode-locked discrete mode diode lasers. These edge-emitting devices incorporate 
a spatially varying refractive index profile which is designed according to the spectral output de- 
sired of the laser. We first demonstrate a device which supports two primary modes simultaneously 
with high spectral purity. In this case sinusoidal modulation of the optical intensity at terahertz 
frequencies can be obtained. Cross saturation of the material gain in quantum well lasers prevents 
simultaneous lasing of two modes with spacings in the millimeter-wave region. We show finally that 
by mode-locking of devices that are designed to support a minimal set of four primary modes, we 
obtain a sinusoidal modulation of the optical intensity in this frequency region. 



I. INTRODUCTION 

Synthesis of terahertz (THz) and millimeter-wave 
(mm Wave) frequencies using optoelectronic means is cur- 
rently the subject of significant interest. Generation of 
these frequencies in the optical domain can be much less 
complex than with electrical schemes and also enables 
signal distribution with the large bandwidth and low 
losses of optical fibers. Interesting applications are ex- 
pected in diverse areas including high-bandwidth wireless 
communication systems, THz imaging and spectroscopy, 
and radio-astronomy P, . 

The most direct route to optical synthesis of these fre- 
quencies is simply to mix the outputs of two detuned 
single mode lasers [11 . This technique is simple to im- 
plement and provides very wide tunability. However, be- 
cause the outputs of the two sources are not correlated, 
broadening of the linewidth of the generated difference 
frequency is observed. For many applications, an active 
feedback mechanism such an optical phase locked loop is 
necessary to improve the linewidth [j, 01 ■ 

An alternative for many applications is to electroni- 
cally modulate the output of a continuous wave single 
mode laser 0] . Electronics based techniques can com- 
bine frequency tunability and narrow linewidth but it is 
necessary to filter the modes of interest from the inten- 
sity modulated spectrum. In addition, higher frequency 
electronic modulation schemes are correspondingly more 
complex and expensive to implement. 

Two-mode diode lasers represent an interesting alter- 
native approach to the problem on account of their com- 
pactness and simplicity. While monolithic two-mode de- 
vices do not generally lead to continuous and wide tun- 
ability of the difference frequency, there may be signifi- 
cant advantages where the frequency of operation is fixed 
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for a particular application. Because the modes share a 
common cavity, the degradation of the linewidth of the 
beatnote due to environmental and other factors is ex- 
pected to be less than in the case of two independent 
single mode lasers Q. Many examples of semiconduc- 
tor lasers that support two modes with a large frequency 
spacing have been described in the literature. These in- 
clude external cavity lasers [9], and devices based on dis- 
tributed Bragg reflections |10l - [l3 | and compound cavity 
designs [Tsl, [l^. 

Here we demonstrate a powerful approach to mode se- 
lection in a Fabry-Perot semiconductor laser that enables 
optical synthesis of frequencies ranging from THz to the 
mm Wave region. Depending on the frequency required, 
two-mode or mode-locked diode lasers that support a 
minimal set of four pre-selected modes can be used. Our 
approach relies on a distributed reflection mechanism 
that is designed as a perturbation of the Fabry-Perot 
mode spectrum. The cavity geometry is derived directly 
from an inverse problem solution that has enabled us to 
design single-mode [11] and two-mode lasers with ex- 
cellent spectral purity. Single mode devices of this kind 
were also shown to have very narrow linewidth emission 
under normal drive conditions [l7| . 

This paper is organised as follows: In section II, we 
demonstrate that two-mode quantum-well lasers provide 
a direct route to optical synthesis of THz frequencies. 
The stability properties of two-mode devices as a func- 
tion of the primary mode spacing are also described. It 
is shown that for primary mode spacings in the mm Wave 
region, cross saturation of the material gain leads to 
wavelength bistability of the primary modes. In section 
HI, we describe how mode-locking of a device designed 
to support a minimal set of four primary modes can be 
used to overcome the problem of wavelength bistability. 
Finally, prospects for obtaining narrow linewidth THz 
and mm Wave signals using discrete mode devices are dis- 
cussed. 
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II. STABILITY PROPERTIES OF TWO-MODE 
DEVICES AS A FUNCTION OF PRIMARY 
MODE SPACING 

We have shown that a set of self-consistent equations 
for the lasing modes can be found by making an expan- 
sion about the cavity resonance condition in a Fabry- 
Perot laser [l^. The effective index along the Fabry- 
Perot cavity is perturbed by N additional features that 
are described by an index step An. A first order scatter- 
ing approximation allows Fourier analysis to be used to 
make a direct connection between the index profile in real 
space and the threshold gain modulation in wavenumber 
space. Here we will present measurements of fabricated 
devices where the index step features are slotted regions 
etched into the ridge waveguide of the laser. All of the de- 
vices considered are multi-quantum well InP/InGaAlAs 
devices with peak emission ranging from 1.3 /im to 1.5 
fim. 

Measurements of a two-mode ridge waveguide Fabry- 
Perot laser with a peak emission near 1.3 /im are shown 
in Fig. [1] (a). Details of the inverse problem solution 
and of the design of two-mode diode lasers can be found 
in [l^. In this case the device length is 350 /im, and we 
have chosen a primary mode spacing of four modes which 
determines a frequency separation of the primary modes 
of 480 GHz. Generally, one finds that as the current is 
increased from below threshold the primary mode on the 
short wavelength side reaches threshold first and as the 
current is increased further the peak power shifts across 
the primary mode spacing to the long wavelength side. 
The spectrum shown in Fig. [T] (a) was obtained at a 
device current of 43.5 mA. At the value of the device 
current shown the time averaged optical power in each 
of the primary modes is approximately equal. We refer 
to this point in the light-current relationship as the two- 
color point. One can see that excellent mode selectivity 
is achieved in this case with the primary mode intensity 
exceeding the background modes by a factor of 40 dB. 

The spectral purity observed in Fig. [T] (a) is limited 
in this case by power transfer to four-wave mixing side- 
bands. The presence of these sidebands is significant as 
they indicate simultaneous lasing of the primary modes. 
This is confirmed by the intensity autocorrelation mea- 
surement shown in the inset of Fig. [T] (a), which confirms 
a sinusoidal modulation of the intensity output at the 
primary mode frequency separation of 480 GHz. Power 
spectra of one of the primary modes and of the total in- 
tensity are shown in Fig. [T] (b). One can see that there 
is significant low-frequency antiphase noise in the spec- 
trum of the primary mode, which can be attributed to 
the global coupling of modes by the carrier density in the 
device. 

These characteristic features of the observed power 
spectra can be reproduced on the basis of a rate-equation 
model for the dynamics of the electric fields of each mode 
El and E2 and for the carrier density N. To account for 
the spontaneous emission noise we add Gaussian white 
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FIG. 1: (a) Optical spectrum of the device at a device current 
of 43.5 mA (two-color point). Inset: Background free inten- 
sity autocorrelation measurement showing mode beating at 
480 GHz. (b) Power spectra of one of the primary modes, ui, 
and of the total intensity, T, measured at the two-color point. 
The dashed lines are theoretical calculations obtained using 
equation \5\ 

noise terms rji (t) and 772 (t) to the dynamical equations for 
the fields and arrive at a system of stochastic differential 
equations of the form 

Ei = ^{l + ia){Ngi-j)Ei+r,i (1) 
E2^U^ + ^a)iNg2-J)E2+m (2) 

N = Jo---NY,9^m\E^\^ (3) 

in 

where the nonlinear modal gain is given by 

ff™=5L°^(l + eE„/3mnl^«P)"'. (4) 

In the above equations, g'^m = 1 is the linear modal gain. 
Jo is the pump current, 7 is the cavity decay rate, and 
Tg is the carrier lifetime. The phase-amplitude coupling 
is given by a. The parameters e/3„m determine the cross 
and self saturation of the gain. A stable two-mode so- 
lution in the free running laser requires /3mn < Pmm- 
The Gaussian white noise terms rim (t) fulfill the relation 
{rjm{t)rjn{t')) — D6{t — t')6mn- A similar model has been 
successfully used to describe the dynamical features of 
two-mode lasers with optical injection [l^ . 

Using the standard theory of multivariate Ornstein- 
Uhlenbeck processes [l^ we obtain the spectrum matrix 
in the small noise limit at the two colour point by using 
the general formula 

S{uj) = ^(J + Iiuj)''^D{J^ ~nu;)-\ (5) 

ZTT 

Here J is the Jacobian matrix of the system [T] - [3] evalu- 
ated at the two-color point and D = diag(Z?, D, 0) is the 
noise coupling matrix. The resulting noise spectra for 
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the individual modes and for the total output are shown 
as dashed lines in Fig. [T] (b). One can see that we ob- 
tain good qualitative agreement with the corresponding 
experimental data, including the characteristic feature of 
the antiphase noise peak at low frequencies of the indi- 
vidual modes and its absence in the total output. The 
analytical treatment reveals that the origin of the an- 
tiphase noise in the dynamics of the individual modes can 
be traced to one negative real eigenvalue of the Jacobian, 
whose eigenvector corresponds to fluctuations around the 
two-color point with a fixed total output and a fixed car- 
rier density. On the other hand, the relaxation oscillation 
peak at around 5 GHz, which is visible in the spectra for 
the total and the individual modes, is related to a pair 
of complex conjugated eigenvalues of the Jacobian that 
describe the exchange of energy between the carriers and 
the fields. For these calculations the value of Jo was fixed 
at 75% above its threshold value, Jthr- Other parameters 
were (7T,) = 800, e = 0.01, and P = Pmn/I^ram = 2/3. 

Time traces of the total intensity and of one of the pri- 
mary modes that illustrate the enhanced antiphase noise 
are shown in Fig. [5] (a). On this scale the total intensity 
is almost constant, while large fluctuations are visible 
in the intensity of the primary mode. These intensity 
fluctuations in the individual modal intensities will lead 
to increased intensity noise in a generated THz signal. 
However, we would like to point out that because the in- 
stantaneous frequency of each mode is determined by the 
dynamics of the carrier density, which in turn depends on 
the dynamics of the total intensity, antiphase noise does 
not lead to an additional increase in the phase noise of 
the THz signal. 

We have found that the level of the intensity noise 
is strongly dependent on how far the two-color point is 
located above the device threshold. By varying the sub- 
strate temperature, the position of the gain peak can be 
moved and the two-color point can be brought further 
from threshold. The normalised standard deviation of 
the intensity noise is plotted in Fig. [l](b). One can see 
that the increased optical power leads to a reduction in 
the level of intensity noise to less than 5%. A device simi- 
lar to that shown in Fig. [T] driven at five times threshold 
was used to generate THz radiation in a photomixing 
setup (21j . 

Because of the increased efficiency of photomixing and 
applications in high bandwidth wireless communications, 
the generation of beat notes in the mm Wave region 
around 100 GHz is of signiflcant interest. However, we 
have found that the simple two-mode device cannot be 
used in this region because of the phenomenon of wave- 
length bistability. Wavelength bistability is expected to 
appear in quantum-well semiconductor lasers once the 
cross-saturation of the gain exceeds the self-saturation 
[l6j . We have found that this boundary is located in the 
region of 300 GHz. For primary mode spacings less than 
this, a discontinuous switching of the intensity from one 
primary mode to the other is observed over a very small 
current range. 
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FIG. 2: (a) Intensity time traces of one of the primary modes, 
v-i, and of the total intensity, T, measured at the two-color 
point, (b) Variation of the normalized standard deviation of 
the intensity noise in one of the primary modes measured at 
the two-color point. Inset: Variation of the device current 
measured at threshold, I^jj, and at the two-color point, Itc, as 
the substrate temperature is varied. 



Optical spectra that illustrate this switching are shown 
in Fig. 131 These are experimental measurements from a 
two-mode device with a peak emission near 1490 nm. The 
primary mode spacing in this case is 250 GHz and one 
can see that a 20 dB change in intensity is observed over 
approximately 1 mA variation in the device current. In 
fact, these quasi-single mode states coexist over a small 
hysteresis region, within which noise driven switching of 
the intensity between the modes can be observed. This 
so-called mode-hopping behaviour is associated with 
a characteristic dwell time, which is the average time 
measured between switching events. By varying the sub- 
strate temperature as before we were able to move the 
bistability region further from threshold. We found that 
the corresponding dwell time varied by over four orders 
of magnitude as shown in Fig. [3] (c). The dashed line 
in Fig. [3] (c) shows a theoretical calculation of the dwell 
time, which is obtained by approximating the system [T]- 
[3]by a one dimensional stochastic differential equation for 

the dynamical variable 4> = arctan / |ii^2|^^ , simi- 

lar to the approach which was taken in [23|. This leads 
to an Arrhenius formula for the dwell time Td given by 



^/2^T 



exp 



1 p2 

2^ 



(6) 



where P = P^-i5/Dts with P = \ [j^ - l) and 5 = 

ie (/3 — 1). For this calculation, we used the same set 
of parameters as in Fig. [T] (b) apart from the ratio of 
cross to self saturation which is now given by /3 = 1.1. 
The noise strength used for fltting experimental data was 
Z? = 9.6 X 10^5 s-2. 
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FIG. 3: Optical spectra of the wavelength bistable device 
before (a) and after (b) the bistability region, (c) Variation 
of the dwell time in the bistability region as the device current 
at the bistability point is varied. The dashed line is obtained 
from the escape time analysis of a reduced one-dimensional 
model [Eqn. [S]. 



III. MILLIMETER- WAVE FREQUENCY 
SYNTHESIS WITH PASSIVELY MODE-LOCKED 
DISCRETE MODE LASERS 
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FIG. 4: (a) Feature density function (solid line). Dashed lines 
indicate intervals where the Fourier transform of the spectral 
filtering function chosen is negative. Inset: Calculation of 
the threshold gain of modes for the laser cavity schematically 
pictured in the lower panel of the figure. Lower panel: Laser 
cavity schematic indicating the locations of the additional fea- 
tures. The device is high-reflection (HR) coated and includes 
a saturable absorber section as indicated, (b) Optical spec- 
trum of the device at a device current of 110 mA and a reverse 
bias voltage of -2.8 V. Inset: Background free intensity auto- 
correlation measurement showing mode beating at 160 GHz. 



In order to avoid the problems associated vi^ith wave- 
length bistability in two-mode quantum-well devices, we 
have designed passively mode-locked devices that sup- 
port a minimal set of four primary modes. These devices 
include a saturable absorber section placed adjacent to 
one of the cavity mirrors. In the device we consider here, 
four primary lasing modes are chosen with a spacing of 
two fundamental modes leading to mode-locking at the 
first harmonic of the cavity. 

The cavity design in this case represents a straightfor- 
ward extension to four modes of the approach to design- 
ing a two-mode device described in JJi] . Fig. |3] (a) shows 
the feature density function which must be sampled in or- 
der to define the effective index profile. A schematic pic- 
ture of the device, high-reflection coated as indicated, is 
shown in the lower panel of Fig. 2] (a) . The device length 
is 545 /im, and iV = 42 etched features are introduced. 
The calculated form of the threshold gain spectrum is 
shown in the inset of Fig. U (a) . 

Experimental measurements of a ridge waveguide 
Fabry-Perot laser fabricated to the design depicted are 
shown in Fig. 2] (b). This device has a peak emission 
near 1.5 /xm and a saturable absorber section of length 
30 /im was placed adjacent to the high reflectivity mir- 
ror. Fig. |l](b) displays the optical spectrum for a reverse 
bias voltage of -2.8 V applied to the saturable absorber 
section. The four selected modes are in excess of 20 dB 
stronger than all other modes and have a frequency sepa- 
ration of 160 GHz. Because most of the laser light is con- 
centrated in the two central primary modes, the resulting 
time evolution of the power is essentially a beating be- 



tween these two modes. The inset of Fig. |4] (b) displays 
the intensity autocorrelation which confirms a sinusoidal 
intensity modulation at the mode-locking frequency. 

We note that while examples of mode-locked diode 
lasers with a sinusoidal intensity modulation have been 
presented in the literature [l^l, the sinusoidal modula- 
tion can be attributed to the fact that the primary mode 
spacing is comparable to the gain bandwidth in these 
devices. For narrower mode spacings a mode locked de- 
vice will generally oscillate on many modes generating 
pulsed output. External filtering of pairs of modes from 
the mode-locked spectrum is then necessary to obtain si- 
nusoidal output [23|. Our approach to this problem can 
be regarded as integrating the functions of mode-locking 
and spectral filtering in a single device. 

These results demonstrate that discrete mode devices 
provide a convenient means to realize THz and mm Wave 
frequencies in the optical domain. Microwave frequency 
synthesis will also be possible in longer devices with a cor- 
respondingly narrower mode spacing. We are currently 
investigating the dependence of the linewidth of the beat- 
note on the drive parameters. Preliminary measurements 
indicate that in the mode-locked device, this linewidth is 
in the region of 5 MHz. It should be noted however that 
because the mode selection mechanism in our device acts 
to predetermine the carrier wave frequency in the mode- 
locked state, we expect that the characteristics of this 
state will also depend on where the central mode fre- 
quency is placed with respect to the material gain and 
loss dispersion. 
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While mode-locking was achieved at 160 GHz in the 
current device, results suggest that quantum-well mate- 
rial can be used for mode-locking frequencies as large as 1 
THz [23|. However, the use of self-assembled gain mate- 
rials may provide significant advantages such as reduced 
amplified spontaneous emission noise, faster carrier re- 
covery and limited alpha factors that can further improve 
the stability of the mode-locked state. The possibilities 
suggested by techniques such as hybrid mode-locking are 
also interesting in this regard. Significantly, single-mode 
discrete mode devices have been shown to have optical 
linewidths as small as 100 kHz under normal drive condi- 
tions [13] ■ This suggests that a mode beating linewidth 
in the sub-MHz range should be achievable with discrete 
mode devices. 



IV. CONCLUSION 

We have demonstrated that discrete mode diode lasers 
provide a direct route to optical synthesis of THz 
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